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ABSTRACT 

The formation of planetesimals requires the growth of dust particles through collisions. Micron- sized particles must grow by many 
orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the 
interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) 
studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a 
previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring 
critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression 
strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and 
conclude that the new code is capable of studying the properties of small aggregates. 
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It is commonly accepted that planets form in protoplanetary 
disks. Giant planets are believed to be produced via one of two 
com peting scen arios. The first one is a gravitational instabil- 
ity dBossl l2007). Clumps of gas (intermixed with dust) may be- 
come gravitationally bound and contract forming giant planets. 
This model is hotly debated (i e JPickett et alJ[200a iBoley et alJ 
2006. 12007b iMaver et al.ll2002l I2004 l200ir The second, com- 
monly accepted sce nario is known as the core accretion model 
( Poll ack et al.ll996h . In this mechanism the solid core of a planet 
forms first and gas is accreted onto that core later. The core itself 
is believed to form by collisional accumulation of planetesimals. 
Terrestial planets, due to their similarity to the giant planets' 
cores, must form out of planetesimals with help of gravitational 
interactions. The formation of planetesimal s is also a subject 
of deb ate. The scenario initially proposed by lGoldreich & Ward! 
(I1973I) assumed a laminar structure of the disk. When the dust 
sublayer reaches a critical thickness, a gravitational instabil- 
ity may form planetesimals from the dense and gravitationally 
bound concentrations of dust. However, this mechanism requires 
an extremely laminar nebula. The shearing motion of the dust 
layer causes a Kelvin-Helmholtz instability and therefore limits 
the thickness of the sublayer. The gravitational instability of the 
dust layer is then prevented . Th is was proved by jWeidenschilling 
dl984l) . ICuzzi etal.ldl993l) and lJohansen et all d2006al) . 

The core accretion model requires to grow micron sized 
dust grains into km sized planetesimals. These are over 27 or- 
ders of magnitude in mass. Radial drift alone may move dust 
particles inwards onto the central star before they reach big- 
ger size. The drift can m ove meter sized particl es all the way 
in within only 100 orbits dWeidenschilUngl fl 9771) . Relative ve- 
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locities of large bo dies set another barrier against accumulation 
of large particles ( Weidensc hilling & Cuzzilll993t I Wurm et al" 



2005a MOrmel & Cuzzill2007l) . A recent study bv lJohansen et af 
(l2006bl) showed however, that planetesimals might be produced 
by a gravitational collapse in the presence of turbulence. Meter 
sized boulders concentrated in high pressure turbulent eddies in 
the disk may form gravitationally bound clumps. Even if this 
process really can be made to work, dust particles must first grow 
over 18 orders of magnitude in mass to be able to form that dense 
accumulations. This growth must happen through the collisional 
sticking of dust particles. 

The initially small, micron sized, dust grains (further referred 
to also as monomers) c ollide and stick to e ach other due to at- 
tractive surface forces d Johnson et al.1 [19711) . The fine dust is 
very well coupled to the gas, meaning that the particles move 
colle ctively. The re lative velocities, due to the Brownian mo- 
tion (lBrown| [l828). are very small. These conditions lead to 
a quasi-monodisperse growth that preferentially collides parti- 
cles of similar size. The aggregates (further referred to also 
as agglomerates or particles) formed in this case are frac- 
tal dKempf et al.ll 1 9991: iBlum et alj|200rl iKrause & Blumll2004t 
iPaszun & Dominikll2006T) ~ 

The Brownian growth produces aggregat es with very open 
structure and the fractal dimension Df - 1 .5 (|Blum et al.l l2000; 
IKrause & Blumll2004 IPaszun & Dominikl 12006). However, the 
gas density influences the growth. When the collisions are no 
longer ballistic, due to high gas density and thus short stopping 
length, the fractal dimension decreases and i n the case of a very 
high g as density it may even approach unity (IPaszun & Dominikl 
2006). This process may play a role in the inner most regions of 
protoplanetary disks. 

The subsequent growth of the dust aggregates eventually 
leads to decoupling of the dust component from the gas. The 
particles st art to be affected by turbulence, se dimentation and 
radial drift (IWeidenschillingl 1 977L [l980l 1 1 9841) . The relative ve- 
locities thus increase and collisions occur preferentially between 



2 



D. Paszun & C. Dominik: Properties of porous aggregates 



agglomerates of different sizes dWeidenschilling & Cuzzilll993l) . 
When the collision energy becomes higher than some threshold 
restructuring energy, aggregates are compacted and their frac- 
tal dimension appro aches Df = 3 (Dominik & Tie lensl 1 1 997t 
Blum & Wurm|p2000l) . Further growth and compaction causes 
the particles to decouple from the g as more effectively a nd the 
relati ve velocities increas e further. IDominik & Tielensl d 19971) 
and lBlum & Wurml d2000l) showed that the outcome of collisions 
can be categorized in terms of collision energy. Perfect sticking 
without restructuring occurs when the collision energy is lower 
than the threshold rolling energy, which is the energy needed to 
roll a single particle over an angle of n/2. When the collision 
energy reaches this limit, the aggregates start to compress upon 
impact. The maximum compaction occurs when the collision en- 
ergy equals the rolling energy times the number of contacts in 
the aggregates. The particles start to lose monomers, when the 
impact energy reaches the number of contacts in the aggregate 
times the energy needed to break one contact. Catastrophic dis- 
ruption occurs when the energy reaches several times the num- 
ber of contacts times the breaking energy. Therefore, as particles 
grow and decouple from the gas, eventually relative velocities 
are re ached that will lead t o shattering of the colliding aggre- 
gates dOrmel & Cuzzill2007l) . 

This general picture has several missing elements. Although 
numerous experiments have been performed i n different size 
and velocity regimes ( Krause & BluirJ l2004t iBlum & Wurml 
200$ IWurm et all l2005albr for a review see Blum & Wurm, 
ARAA,2008 in press), our understanding of the processes in- 
volved in the growth of aggregates is still incomplete. The re- 
structuring mechanism for instance is understood only qualita- 
tively. The degree of the collisional compaction still remains a 
mystery. This, however, is crucial for the growth of the meter 
sized aggregates. The sticking efficiency as a function of the par- 
ticle density along with the density evolution must be determined 
in order to understand quantitatively the growth of meter sized 
boulders. 

The low strength of the aggregates due to the fractal struc- 
ture leads to fragmentation in the case of fast impacts. The dis- 
tribution of fragment sizes must be also determined as a func- 
tion of the collision energy. This makes it possible to keep 
track of the realistic size distribution in the disk. Small par- 
ticles, if not replenis hed in the disk, are very qui ckly swept 
up by larger grains (Dulle mond & Dominik! 120051) . Thus the 
small particles should be supplied to the disk by some pro- 
cesses. Dominik & Dull emondl (|2()08) show that the infall of 
grains from the parent cloud is rather unrealistic and requires 
fine tuning of parameters to reproduce observational data. Thus 
collisional fragmentation is the most likely mechanism that can 
explain the population of small grains in protoplanetary disks. 
However, long before fractal aggregates approach velocities high 
enough to cause fragmentation, they undergo collisional com - 
pression dBlum & Wurdl2000tlWeidenschilling & Cuzzilll993l) . 
Thus the only scenario leading to fragmentation of fractals is 
collision with much larger non-fractal particle. 

The current understanding of processes like sticking, bounc- 
ing and fragmentation of aggregates is poor. The understanding 
of this processes on micro scales may fill these gaps and allow 
extrapolations to the larger sizes. This may be resolved by two 
approaches. Experiments may be performed in a laboratory to 
provide useful data. However, this way is available only for ag- 
gregates below centimeter in size. Larger particles are currently 
inaccessible for experiments. The second, theoretical approach, 
understanding the material properties of porous matter reach ag- 
gregates of meter size and beyond. Thus it would be ideal to 



provide theoretical predictions for centimeter sized and larger 

ag gregates 

ISironol d2004l) developed a model capable to simulate big 
- meter sized and larger aggregates, using Smoothed Particles 
Hydrodynamics (SPH). In this case one particle in the model is 
an aggregate characterized by material properties like, compres- 
sive strength, tensile strength, density and sound speed. In order 
to obtain some of these properties he fitted power-law functions 

to experimental data of compression and t ensile strength . 

This method was then also used by Scha fer et al.l d2007l) . 
Further studies require, however, a more precise determination 
of the material properties of porous bodies. These may be ob- 
tained in laboratory experiments or in computer simulations as 
presented in this work. 

More real istic Solar Nebula dust analogs were used in ex- 
periments by IBlum et all d2006l) . They investigated aggregates 
made of micron sized dust particles. Moreover different aggre- 
gates used in these experiments consisted of spherical or irreg- 
ularly shaped monomers. The compression and tensile strength 
curves for this dust cakes were determined. 

Although experiments give a quantitative description of pro- 
cesses in cm sized aggregates, it is difficult to access for small, 
a few microns size aggregates. Aggregate dynamics models 
dDominik & Tielensl \\Wft 



Dominik & NiiboL 



dynamic 
<f l2002l) 



are ideal 

for simulations of these small scale structures. This method spa- 
tially resolves single monomers, which is certainly needed to un- 
derstand physics of bigger - meter sized bodies. Until now the 
main drawback of this aggregate dynamics model was missing 
quantitative agreement with experiments, even t hough the qual- 
itative agreement has already been established (Blu m & Wurml 
2000). Thus a calibration of this model is required. 

Another aggrega te dynamics model was recently presented 
bv lWada et al.ld2007l) . They made use of potential energies in or- 



der to derive forces acting between grains in different degrees of 
freedom. In this case they present just a 2D case, but their results 
are in agreement with findings of N-body model presented ear- 
lier bv IDominik & Tielensl d 19971) . Wada et al. (ApJ,2007; sub- 
mitted) shows the first results of compression and disruption of 
3 - dimensional aggregates in head - on collisions. Although 
the model is qualitatively in agreem ent with experiments, as is 
studied by Dominik & Tielens ( 1997), the quantitative mismatch 
is still present. IWada etal.ld200l " do not show any solution or 
workaround to the quantitative disagreement between theory and 
laboratory experiments. The work presented here addresses this 
issue and provides mechanisms that fit our model to the empiri- 
cal data. 

In this paper we present t he calibration of the aggre gates 
dynamics model develope d bv IDominik & TieTe ns (1997) and 
Dominik & Niibold (2002). We fit the code using experimental 
data and further test it. Modifications that are implemented in 
the code are presented together with a few possible application 
of the model. 

2. The model 

In order to study the agglomeration mechanisms involved in the 
growth of planetesimals, we use the SAND code (Soft Aggregat e 
Numerical Dynamics) developed by Dominik & Niibold (2002). 
It is a N-body model of a system of spherical particles interact- 
ing via surface forces. 

Two monomers feel the attractive force only when they are 
in contact. The surfaces of the particles deform and form a con- 
tact area. When the particles are pulled outwards, increasing 
the relative distance, the area decreases and the monomers are 
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pulled back to the equilibrium position by the surface forces. 
The compression of the system on the other hand leads to in- 
crease of the contact area and repulsive force pushing the parti- 
cles apart . A detailed descript ion of the surface forces was pro- 
vided by Johnson e t alj dl97ll) (further referred to as JKR). The 
influence of adhesi on forces on the conta ct between particles 
was also studied by Deri agu in et al.l d 19751). The model is abl e 
to treat long range magnetic forces (Dominik & Nubold 2002), 
however these are not the subject of our current study. 

There are two main processes that govern the events dur- 
ing a collision. The first one is breaking a contact between two 
monomers. This dissipates part of the energy and weakens or de- 
stroys aggregates participating in the collision. The second pro- 
cess is a rolling motion of a monomer over another grain. This 
also dissipates energy and causes a restructuring of an aggre- 
gate. This restructuring may be attributed to both compression, 
which results in strengthening the aggregate, or decompression, 
i.e. weakening the aggregate's structure. 

The JKR theory predicts a critical force that is needed to 
separate two particles. This p rediction was test ed for the case of 
micron size Silica spheres by lHeim et alJ (1 19991) . Two monomers 
were pulled off each other using an Atomic Force Microscope. 
The measured pull-off force was in agreement with the force pre- 
dicted by the JKR theory. The pull-off force is 

F c = 3nyR, (1) 

where y is a surface energy and R = (j^ + -^j is a reduced ra- 
dius of the spheres in contact. Thus the results of the experiments 
confirm the theoretical predictions. 

A similar experiment was performed in order to determine 
the horizontal force s acting on the particles in c ontact (first stud- 
ied theoretically by Do minik & Tielensl dl997l) ). The horizontal 
displacement of the contact zone causes a torque acting against 
the displacement. The resulting rolling frictio n was measured i n 
laboratory experiments, again using an AFM (iHeim et al.fl 999). 

2. 1 . Rolling friction 

Rolling friction is one of the most importan t energy dissipation 
chann els in the restructuring of aggregates ( Do minik & Tielensl 
I1997I) . It is thus of a great im portance to treat it in a corre ct way. 
The theoretical derivation by Domini k & Tielensl dl997l) shows 
that the energy associated with initiating rolling is expressed as 

e ro0 = 6nyf ctiv (2) 

where ^ crlt is a critical displacement at which the inelastic behav- 
ior occurs and energy is dissipated. £ CI j t was initially assumed to 
be of the order of inter - atomic distances dDominik & Tielensl 
I1997I) . The rolling motion causes shift of the contact area, im- 
plying that at one end the contact has to be broken in order to 
form it at the other side of the contact area. The experimentally 
determined friction showed that the critical displacement must 
be larger. The determined value was £ cr ; t = 3.2 nm meaning that 
the dis placement is rather of the order of ten inter atomic dis- 
tances. Dominik & Nubold (2002) took that already into account 
and used in their model the value of f cr j t = 10 A. This however 
was applied to diffe rent material th an the one investigated in the 
lab. We followed IHeim et all d 19991) and applied larger value of 
ftrit = 20 A, which is approximately 10 inter atomic distances. 



2.2. The pull-off force and the critical velocity for sticking 

At low velocities the dominating mechanism that dissipates 
energy is rolling. At higher impa ct speeds, however, other 
channels become more important. IChokshi et al.l (1993) and 
Dom inik & Tielensl (1 19971) calculated how much of the initial 
energy is dissipated in the collision between two particles. This 
gives the maximum energy at which the particles can stick in a 
head-on collision. The critical velocity is given by 

v -t = l-07 £ . 1/3j?s/6pl/2 , (3) 

where y, R and p are surface energy, reduced radius and mass 
density, respectively. E*~ l = -gp- + i^ 2 - is a reduced elastic- 
ity modulus with v, and Ej being Poisson's ratio and Young's 
modulus, respectively, of grain i. When eq. [3] is applied to a 
R\ = 0.6 jum silica grain, impacting flat silica surface, it gives 
the critical v elocity v cr it = 0-1 8 ms _1 . This was again tested in 
experiments. Popp e et ail d2000h showed that such a particle can 
stick to the surface at significantly higher velocities, of the or- 
der of 1 ms _I . The measured velocities were 1.2 + 0.1 ms -1 for 
0.6 jjm grains and 1 .9 + 0.4 ms _1 for 0.25 fim grains. They used 
slightly different definition of the critical velocity. The critical 
velocity was defined as the velocity at which the probabilities of 
sticking and bouncing are equal 50%. They measured the stick- 
ing velocity for different materials, shapes and sizes of particles. 
The resulting critical velocity was shown to depend on size of 
the monomer as a power law with index of about 0.53. Although 
the theoretically derived slope R~ 5 t 6 of the power law is differ- 
ent from the empirical data, they are still consistent within error 
bars. Moreover the power law was fitted to only two data points. 

We will assume that the theoretical dependence on the ra- 
dius is correct and consistent with the experiment. The discrep- 
ancy between the absolute values of the critical velocity however 
have to be revised, doing so is absolutely essential for meaning- 
ful results. 

In order to better understand the mechanisms inv olved in the 
stickin g of the monomers we refer to experiments by Popp e^eTail 
(2000) and iDahnekd d 19751) . They investigated the bouncing of 
micron sized spheres off a flat silica surface. In the former exper- 
iment the impacting particle was made of silica, while the latter 
on e used soft polyst yrene grains. 

iDahnekj (1 19751) determined experimentally the critical ve- 
locity to be about 1 ms _1 . The ratio of the rebound velocity to 
the incident speed (coefficient of restitution) decreases in this 
case with decr easing impact v elocity. At very large velocities of 
about 20 ms _I Dahneke d 19751) noticed a decrease of the resti- 
tution coefficient with increasing impact speed. He related that 
behavior to plastic deform ation. 

In the experiment by Poppe etail (120001) the first positive 
value of the coefficient of restitution indicates a critical veloc- 
ity of the order of ~ 1 ms" 1 . 

In order to scale the model to the experimental data we in- 
troduce a mechanism that dissipates part of the initial energy at 
the moment of first contact. The plastic deformation of small as- 
perities on the surface of the grain is an easy way to dissipate the 
energy and increase the critical velocity. Below we will estimate 
the central p ressure in order to compare that to the yield stren gth 
of 104 MPa (ICallisterll2000l) for siliciPl lChokshi et all dl993l) on 
the other hand argue that the yield strength for very small bodies 



1 Since the yield strength of brittle material like silica is not defined, 
flexular strength is given. 
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is of the order of 0.2 times the Young's modulus of the material. 
This corresponds to 1 1 GPa for silica. 

The maximum pressure in the contact area occurs in the cen- 
ter of the contact. The rad ial distribution of pressure in the con- 
tact zone is the following ( Domini k & Tielensll 19961) 



P (r, a) = 6 ^ * (1 - m> - 2 *L( * Pd - (^V * (4) 
ao in this equation is the equilibrium contact radius given as 

p2\l/3 



l9nyR' 



(5) 



Now we need to estimate the maximum contact radius that is 
reached during the collision. The radius of the contact area in- 
creases with increasing normal load. The force applied to the 
sphere during the collision is F — d(mv)/dt, where d(mv) is the 
momentum of the impacting particle and dt w 1CT 9 s is the col- 
lision time. With this force applied, the contact radius is 



a — 



I 3R 

4£* 



(F + dnyR + ^(6nyR) 2 + UnyRF) 



1/3 



(6) 



as derived by lJohnson et alJ (Il97lh . Thus in the case of the ve- 
locity v = 1 m/s, the central pressure is of the order of 1 GPa. 
This pressure however is exactly in between the two values of 
the yield strength specified above, making the plastic deforma- 
tion possible. The lower limit of the yield strength (104 MPa) is 
reached at radii r/a < 0.914, which makes 83% of the contact 

area exposed to potential plastic deforma tion. 

This idea was investigated before bv lTsai et al.l d 19901) . The 
energy consumption in such a deformation can be expressed as 



J asp 



YV. dsp , 



(7) 



where Y is the yield strength of the deforming material and V asp 
is the volume of the asperities that are flattened during the colli- 
sion. To calculate the dissipa ted energy we ne ed to get the vol- 
ume of asperities. We follow ITsai et al.l d 19901) again and calcu- 
late the maximum contact area given by 



' 3R 



[F eq + 6tt 7 R + ^(6nyR) 2 + \2n jRF eq ] 



1/3 



(8) 



(IChokshi et al.ll 19931) and the equivalent impact force given by 
ITsai et all d 19901) as 



F eq =2.5 3/5 4^ 1/5 £* 2/5 . 



(9) 



Now we can get the volume of the deformed material by multi- 
plying the contact area by the volume of a single asperity, and 
the number of asperities per surface area 



V aS p = -T^aspNasp 



(10) 



where r asp is the radius of a single asperity. Here we assume the 
bumps to be hemispheres distributed homogeneously over the 
surface of a particle or a target. In our case we use a parameter 
that describes the total volume plastically deformed per unit of 
the surface area. This way we have only one parameter that de- 
fined how efficiently the energy is dissipated. Note that in reality 
the pressure in the contact area is compressive in the center and 
tensile on the edge. Thus our parameter may be slightly higher 
than what we present. 
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Fig. 1. Critical velocity as a function of a grain size. Squares 
indicate the model without additional dissipation process, dia- 
monds show present model with the energy dissipation by plastic 
deformation of surface rou ghness and triangle s show experimen- 
tally determined values bv lPoppe et alJ |2000) with error bars. 



Fi g. [TJ shows the cr itical velocities determined experimen- 
tally (Poppe et al. 2000) for two different grain sizes. Also the 
critical velocity obtained using our model is plotted showing that 
we successfully fitted our model and we can reproduce the ex- 
perimental results. 

In our model the volume of the asperities deformed upon 
collision is very small. When we assume surface roughness to 
be hemispheres with radii r asp = 1 nm, the fraction of the area 
occupied by the asperities is only about 2%. We can thus still 
say that the molecular size asperities may be responsible for in- 
creasing the critical velocity. Moreover we see in fig. Q] that the 
dependence of the sticking velocity on the grain size is not af- 
fected and the only difference is that the energy regime has been 
shifted to the level measured in experiments. 

In order to avoid confusion, we must stress here again that 
we do not claim that the plastic dissipation takes place during 
the collisions of silica spheres. We simply implement additional 
scaling parameter in order to fit our model to the experimen- 
tal data. The discrepancy between the empirical data and theory 
should be further investigated. We also think that plastic defor- 
mation might need somewhat more attention, because of the very 
high pressure that is present in the center of the contact area and 
relatively low strength of the material considered in this work. 

2.3. Excitation and Cooling 

In parallel to the energy dissipation due to contact breaking, en- 
ergy can be dissipated via other channels in the vertical degree of 
freedom (along the line connecting centers of two grains in con- 
tact). Two grains held together by the surfac e force vibrate rela- 
tive to each other (Dominik & Tielens 1995). Ideally the oscilla- 
tion is frictionless and no dissipation occurs. In reality, however, 
the vibration causes an oscillation in the size of the contact area. 
When decreasing the contact size, part of the energy is dissipated 
by breaking the connections at the edges of the contact area. This 
ultimately leads to a cooling of the aggregate and damps the vi- 
brations. If this mechanism is not taken into account, successive 
slow collisions may heat up the aggregate and ultimately lead to 
"evaporation" of monomers from the aggregate surface. 
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In fact we observed this phenomenon in our simulations of 
linear chains. Successive collisions of grains with an aggregate 
caused an increase of the amplitude of the oscillation and eventu- 
ally lead to breaking of several contacts. All individual collision 
velocities were far below the sticking velocity, proving that this 
mechanism may produce artificial results. 

To resolve this problem we introduced a weak damping 
force, which was intended to slowly dissipate the vibrational en- 
ergy. The force is acting in the vertical direction in respect to the 
contact area. The damping force is expressed as 



iamp 



const V, 



(ii) 



where v- is the vertical component of the relative velocity, and 
const is arbitrarily chosen to damp exponentially 95 percent of 
the vibration energy within ~ 100 oscillation periods. By tuning 
the constant to this value we made sure that the damping force 
influence is not significant within short timescales of a few first 
vibration periods, where the sticking or bouncing event is deter- 
mined. In the first period the dissipation takes place and if the 
particle looses not enough energy it will bounce. If the damping 
force is too strong, it may remove enough energy to allow stick- 
ing. We made sure, that the main energy dissipation process that 
sets the critical velocity is the plastic deformation mechanism 
and the sticking velocity is not affected by the presence of the 
damping force. 

To show the energy leakage due to the damping force we 
calculated the total energy in a vibrating pair of monomers. The 
particles were displaced from the equilibrium position and the 
kinetic and potential energies were calculated. In the case of a 
frictionless oscillation the potential energy is 



Jo 



F(6')d8', 



(12) 



where 6 is a displacement such that the distance between centers 
of the two grains is r\ + - 6. The vertical force 



F(a) 



4F c ||iL 



). 



3/2 



(13) 



depends on the contact size a and equilibrium contact radius «o- 
The integral (fT2l may be then changed into 



£ p = J F{a')da. 



(14) 



To get the energy we solve eq.[14]by changing the variable from 
8 to a using the following relation 



2 R \ U 



3\ao/ 



The potential energy is then given as 

,5 „ / x 7/2 



a l 

4F C — 
R 



2(aY _ 21 aY 
5\ao) 3\a ) 



1 a 
6\flo 



(15) 



(16) 



The kinetic energy is simpler. We just add kinetic energies of all 
monomers together 



~2~ 



(17) 



where m, and v, are mass and absolute velocity of the 2-th grain. 

The total energy E lot in the case of a frictionless oscillation 
is plotted in figj2]a. To better see variations in the kinetic energy 



we shifted it to the level of the potential energy. The shift is 
equal to the potential energy of the system in the equilibrium 

£ p = jr° F(5')dd'. The total energy in this case is conserved and 
the amplitude is constant. When we enable the damping force 
the energy starts to leak. The potential and kinetic energies may 
be calculated using the same formulas. The energy dissipation is 
presented in fig.[2]b. 

The energy loss due to the damping force is very weak and 
removes about 95% of the total energy within approximately 100 
vibration periods. Thus other processes like rolling, and break- 
ing dominate the energy dissipation. In particular the critical ve- 
locity for sticking is not influenced at all by this damping force, 
because the assumed plastic deformation of asperities dissipates 
nearly the entire collision energy within the first vibration pe- 
riod. 



3. Applications 

The presented model is a very good tool that can be used to 
study the aggregation of dust particles. However, the modifi- 
cations that we introduced in the previous section need further 
verification. Although the critical velocity is well fitted to the 
experiments, it is very useful to further test the model against 
laboratory experiments. When the model is successfully verified, 
it can be applied confidently in the study of aggregates dynam- 
ics. Compaction and fragmentation of dust aggregates are pro- 
cesse s that require further investigation ( Dullem ond & Dom inik 
120051) . Our model is perfectly suited to this task. Dynamical pro- 
cesses, involving small sub-mm sized aggregates, can be studied 
in depth using our tool. We can provide detailed understanding 
of micro-physics that governs compaction and fragmentation in 
this scale. Phenomena involving aggregates of mm and larger 
sizes require an entirely different approach. They must be treated 
with different methods i.e. Smoothed Particles Hydrodynamics 
(SPH). However, to do this, material properties, such as com- 
pressive strength and sound speed, are needed. Below we present 
simulations which can be directly compared to experiments and 
provide great test of our model. We simulate compression of a 
dust cake and determine the compressive strength of such cake. 
We also determine the sound speed in such a porous aggregates. 

3.1. Compression 

Blum & Schrapler (2004) measured the compressive strength of 
dust cakes in the laboratory. The sample was prepared by ran- 
dom ballistic deposition (RBD). Single monomers were shot 
from one direction at low velocity (hit-and-stick) and grow the 
dust agglomerate. The resulting "cake" was about 2 cm in di- 
ameter and similar in height. The finished cake was later placed 
between two flat surfaces. The load was applied onto the upper 
surface and caused it to move towards the lower surface com- 
pressing the dust sample. In order to simulate compression of a 
dust cake we need to prepare adequate setup. 

3.1.1. Setup 

In the experiment bv lBlum & Schr apler (2004) the applied pres- 
sure resulted in compression of the dust cake. Measurement of 
the cake volume resulted in determination of a filing factor (f>. 

Our setup was organized in a similar manner to the experi- 
mental one. Our code, however, can handle only spherical par- 
ticles. Thus instead of two planes we used two very big "wall" 
grains, each with diameter of 2- 10~ 2 cm. The sample was shaped 



6 



D. Paszun & C. Dominik: Properties of porous aggregates 



1.624 
1.626 

1.628 



1 .632 



£ -1.634 



1 .636 
1 .638 



a 



<\ ■"■ /, .'■ <\ ■•■ :■ i M 



/ .-■ m ;', / ' r 



ii ■■ : i -■; / , ; M . '1 ;■ l\ 



1 ii ,i i! l. ; i' ii l. : i' ij ,i Ii 

'j 



1 >- : ; ; i f ■ t I 1 ii i. : i 1 ; i. ; ; ' t i. : ':' 



t I j ! I ! i i ! j I ■ | ! i ! 
f : I i ? f . f j ''• .1 -l r : -I t 
'i , :i i: :', '': : > 'i ;1 i': : i : ; ii ,: :i '': : i 




,:' i' l:' < \l \ ,: » : ■ 

! i ! * ! ? ? K n >; i 1; : 

'! ii ••' ,! i 1 i 1 :'l i: S M. 



1 2 3 4 

time [10" 8 s] 



12 3 4 5 

time [10" 8 s] 



Fig. 2. The total energy is redistributed between the kinetic energy (dashed line) and the potential energy (dotted line). At the 
maximum and minimum separation the potential energy is maximum, while in the equilibrium position the kinetic energy is the 
highest. The kinetic energy was shifted down for better overview (see text). No damping case - a, and the energy leak form the 
system due to the damping force - b. 



as a cylinder and the distance between two compressing "wall" 
grains was adjusted to exactly fit the micro dust cake. 

Our dust cake was grown via particle cluster aggregation 
method (PCA). The target aggregate (initially a monomer) was 
randomly oriented before each collision with a grain approach- 
ing at a random impact parameter. Any successful hit resulted 
in perfect sticking without any restructuring. The new target was 
then randomly oriented, and a new grain was shot at it again. The 
resulting aggregate was then shaped into a cylinder by removing 
all grains outside of the desired contour. 

The size of the cake was 13.2jum in diameter and it was 
13.2jum high. When filled by 291 monomers with radius of 
0.6 yum the filling factor of the cake was <p = 0.146. For compar- 
ison the dust cake used in real experiment had the filling factor 
4> = 0.15 + 0.01. The monomer size in that experiment was also 
of similar size ro = 0.75 //m. 

The compression was done by moving one of the large grain- 
surfaces at constant velocity. While it was approaching, the sec- 
ond grain was fixed in its position and was unable to move even 
upon extreme pressures. The sample was thus compressed by 
one surface against the other one. The initial setup is presented 
in fig. [3^. Two large grains on both sides of the aggregate are 
the back wall, fixed plane in the right and the approaching, com- 
pressing plane in the left. The aggregate is placed in between 
and the particles can escape sideways increasing this way the 
diameter of the cake. 

In order to simulate quasi static compression we fixed the 
velocity of the compressing "wall" grain to 0.05 m/s. This is 
over an order of magnitude lower than a critical velocity and 
much lower than the sound speed in this medium (see section 
3.2). Thus the assumption that we are in a quasi static regime is 
reasonable. The dominant acceleration of a single monomer, in 
this case, is due to surface forces. 

For the purpose of this compression simulation, we disable 
the net force acting onto the compressing "wall" particle, but 
we save the record of this force for each time step. Thus the 
approaching surface cannot be stopped and moves with constant 
speed. For each time step the net force is stored and later used to 
determine a pressure. 



At each successive step the dust cake becomes slightly more 
compressed. The degree of compression can be related to the 
force that was needed to get the cake into this state. 

3.1.2. Results 

The small size of the dust cake used in this numerical experi- 
ment causes a low number of contacts of the compressed aggre- 
gate with the approaching surface. Consequently, the net force 
applied to the compressing surface was strongly variable. Every 
new contact formed between the dust cake and the incoming sur- 
face resulted in a sudden decrease of the force. The new con- 
nection is initially stretched and thus the surface is attracted. 
Similarly, oscillations of monomers at the surface cause addi- 
tional variation. This makes it more difficult to uniquely deter- 
mine the compression force. To overcome this problem, for ten 
each successive time steps we choose the one with the highest 
force, because ultimately this is the force required to compress 
the cake. 

The pressure was calculated as P = j, where F is the normal 
load and S is the cross-section area. 

Fig. [3] shows the initial setup and the results of compression 
with increasing pressure. The lowest pressure cannot restructure 
the aggregate. The height of the dust cake is almost unchanged. 
A higher pressure of 2 kPa compresses the cake. The monomers 
in the cake's surface are pushed down into the cake. At a pres- 
sure of 1 • 10 4 Pa the cake is compressed significantly, causing 
a horizontal flow of the particles and thus an increase of the 
cake diameter. The evolution of the dust cake cross-section is 
shown in fig. [4] The relative increase is affected by the size of 
the dust cake. What may be a negligible boundary effect in a 
large macroscopic aggregate, here it has a large impact onto the 
entire dust cake. When the diameter of the cake increases just by 
a diameter of a sin gle monomer, the area inc reases by 40 %. In 
the experiment by Blum & Schrapler (2004) a cake of about 2 
cm diameter was compressed and the final increase of the cross- 
section was measured to be larger by a factor of 1 .6 relative to 
the initial cross-section. In our simulation the area increases al- 
most 4 times, subject to significant uncertainties in the deter- 
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Fig. 4. The ratio of the cross-section of the dust cake to the initial 
cross-section as a function of the applied compressive pressure. 
The data presented here is from a single experiment. Each point 
is plotted with a 40% uncertainty in determination of the cross- 
section (see text). 
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Fig. 5. The volume filling factor vs normal pressure. The 
solid line indicates the re sults of the laboratory experiment by 
Blum & Schrapler (2004). The dashed area is the error to that 
data. Squares show results of our single simulations. Error bars 
are due to difficulty in determination of the volume of the final 
aggregate. 



mination of the dust cake cross-section. The initially cylindrical 
shape changed upon compression into an irregular profile. 

To determine the volume filing factor, we used only the inner 
part of the cake. A cylindrical volume enclosing initially the dust 
aggregate, and used to determine the filling factor, was decreased 
in height only. In this way we reduced the uncertainty that arose 
from the boundary effects. The initially porous aggregate is de- 
formed and can expand. Particles are pushed into the cake, fill- 
ing voids. Thus the filling factor must increase as an effect of 
compression. Fig.|5]shows our results together with the data ob - 
tained in the laboratory experiments (Blum & Schrapler 2004). 
Low pressures are unable to affect the aggregate. However, the 
boundary effects also cause problems. Initially the "wall" parti- 
cle connects to the cake by only one monomer. This causes very 
strong variations in pressure. We show the compression curve 
starting in the point when the "wall" particle has made 5 con- 
tacts with the cake. At this point, however, a few grains are al- 
ready pressed into the dust cake, causing an increase of the vol- 
ume filling factor. Thus our compression curve in fig.|5]is shifted 
upwards. The filling factor is also overestimated by using big 
spherical grains as the compressing surfaces. The height of the 
dust cake we used is calculated as D - 2/? compress , where D is dis- 
tance between centers of the two big grains and /Compress is their 
radius. Thus the volume occupied by the dust cake is actually 



slightly larger. This effect is initially less important, as it con- 
tributes only about 2 % error. However, it gets more relevant at 
larger pressure, where the cross-section of the dust cake is larger 
and the volume smaller. 

Error estimation was done using the central parts of the dust 
cake, where boundary effects are smaller. The filling factor was 
determined by calculating the volume of monomers enclosed in 
a cylindrical volume. The error was then estimated to be the dif- 
ference between the filling factors determined for two different 
cylinders. The first one had a radius of 6 fim, while the second 
one was one monomer radius smaller. 

The volume filing factor is initially constant, until the pres- 
sure reaches the value of about 5 • 10 2 Pa. The filling factor 
increases until it reaches the maximum compression of about 
<f> = 0.35 and remains constant. The most interesting result is the 
resemblance of our findings to the lab expe riments. The onset 
of com pression fits the experimental data of Blum & Schrapler 
(2004I). Our compression curve follows the laboratory data 
tightly ending up at a similar value of the filling factor. The small 
differences between two curves are most likely due to small size 
of our simulated dust cake. 
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3.2. Sound speed 



Table 1. Material properties used in this work for silica. 



One of the very important properties in the porous material is 
the sound speed. It must be lower than the bulk sound speed be- 
cause the mass is being moved by a force acting on a very small 
contact area. The collision of two grains with supersonic veloc- 
ities can result in complete disruption of those bodies. We will 
first derive an analytical formula for the sound speed. For this we 
apply the JKR theory and assume that the signal transported is 
a very small perturbation. The main assumption is that the two 
monomers in contact behave like a perfect spring which is the 
case for small amplitudes. However, the contact forces are asym- 
metric in respect to the equilibrium position. Thus compression 
of two particles results into different forces than stretching it to 
the same displacement. We can therefore expect that large am- 
plitudes lead to modified sound speeds. In the following sections 
we present an analytical approach and later we show our simula- 
tions for both the simplified case of a linear chain of monomers 
and the general case of a non fractal aggregate. 

3.2.1. Analytical solution 

Every two monomers that are in contact are held together b y 
surface forces acting in the contact area dJohnson et aljfl97ll) . 
Mutual attraction inevitably leads to a vibrational spring-like 
motion, which in linear approximation can be written as 



-k(6 - 6), 



(18) 



where 6, do, k and F are displacements, equilibrium displace- 
ment, spring constant and restoring force, respectively. Note that 
the displacement 5 is defined in a way that it increases, when the 
distance between two monomers is decreasing. In our case we 
want to determin e the spring constant, which is necessary to find 
the sound speed. Johns onet al.l (Il97ll) showed that the surface 
force in the contact area is related to the contact area radius a by 



W(±)'-(f v: 

\\aoJ \a 



(19) 



We also know the relation between the displacement and contact 
radius 



6 = 6 I/3 5 c (2a V -4/3a I/z fl 1/z ), 



-2.2 



(20) 



where 6 C = l/2gf^ is a critical displacement for breaking the 
contact. With these equations we can determine the spring con- 
stant and later the sound speed. First we differentiate eq.[T9land 
l20lat a — «o to get 



&F_ 
da 

and 

d&_ 
da 



6F C 
a 

fl() ' 



(21) 



(22) 



where So = is an equilibrium displacement. Now we can write 



dF_ _ d£da _ dF/ddy 1 

d6 da d6 da \ da I 

and substituting eq.|2T|and[22]we get 

dF = ^(9nyR 2 E* 2 ) 1/3 dS. 



(23) 



(24) 



E* [dyn/cm 2 ] 


y [erg/cm 2 ] 


Po [g/cm 3 ] 


2.78 • 10" 


25.0 


2.65 



Since the restoring force is exactly opposite, we may add the 
minus sign here and see that the spring constant k is 

. 6 



;(9nyR i E i " L ) 



2 f *2n1/3 



(25) 



With the spring constant of the oscillating system of two 
spheres in contact we can proceed to the calculation of the sound 
velocity. For a spring the velocity of sound is given by 



Cs = 



(26) 



where L is the total length of the spring (or set of the springs) and 
pi is a mass of the spring per unit length. We can then determine 
the velocity in a linear chain of n grains in contact, each ro in 
radius. The length of such a string of grains is 



L = n(2r - S ) + S - 2r . 



(27) 



Since the spring constant for a set of springs is simply = k/n, 
the sound velocity can be written as 



c s = 



k L L kL 2 /n 



nm n V 4/37ir 3 po ' 
Substituting eqs. [25] and [27] we get 



1\ (9nyR 2 E* 2 )^6/5 

1 - - (2r - So) J — ; — • 

n) V 47T^ P0 /3 



(28) 



(29) 



Now we can try to see how the sound velocity in such system 
depends on the size of a single monomer. 



c s oc ro 



2/3 



.-1/6 



(30) 



This means that to double the speed of sound we need to use 64 
times smaller grains. Therefore, the sound speed in this case is a 
very weak function of monomer size. 

The sound speed of an infinitely long chain of 0.6 yum silica 
grains is c s = 513 m/s. If we apply our findings to a chain of 
50, 0.6 micron sized, silica grains, the sound velocity turns out 
to be c s = 503 m/s. We can compare it now to the previous 
results obtained in research of granular medium composed of 
mm sized and bigger grains. In all our simulations and analytical 
cal culations we u s ed mat erial properties as specified in tab. [TJ 

lHascoet et al.l ( I 19991) developed a model of macroscopic 
grains to study the propagation of sound in a granular chain. 
The centimeter sized grains they use do not interact in this case 
via attractive surface forces. The signal is transported due to 
the Hertzian stress that arises as an effect of overlap of grains 
in contact. They also apply the spring theory but with an arbi- 
trarily chosen spring constant k. For values of k in the range 
between k = 10 6 N/m and k = 10 8 N/m, and mass den- 
sity of 1.9 x 10 3 kg/m 3 , the sound velocity they derive is in a 
rang e of 300 m/s t o 300 m/s. Similar results were obtained 
by iMouraille et al.l d2006l) . A sound speed c s » 200 m/s was 
found for closely packed grains. In this case monomer radius 
was 1 mm, density 2 x 10 3 kg/m 3 , and the spring constant 
k = 10 5 N/m. 
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3.2.2. Numerical experiment 

In order to verify our findings we performed a numerical exper- 
iment. We prepared a linear chain of 50 monomers. The first 
one in line was slightly displaced from it's equilibrium posi- 
tion. When the simulation started the grain started to move in 
order to reach it's equilibrium and the second grain was dis- 
turbed. Such motion was propagating until it finally reached the 
last grain. The total distance traveled by the density wave was 
99(ro - 6o). For 50 silica monomers with radii ro = 0.6 fim 
the travel time took t — 1.16 10~ 7 s, which results in the sound 
speed of c s = 512 m/s. We can now compare the value with the 
theoretically derived sound speed c s = 503 m/s. The difference 
between the theoretical velocity and the one that was obtained 
numerically is only about 3%. Moreover in the simulation the 
first grain is displaced by a finite distance meaning that the as- 
sumption of low displacements may not be entirely correct. The 
displacement of about 0.5 6o is relatively large. We performed 
a series of simulations with stronger perturbations and the data 
obtained shows that the sound velocity increases as the pertur- 
bation strength increases. 

3.2.3. Porous aggregates 

Since a linear chain of monomers is a special case, we will 
discuss now the more general agglomerates of irregular shape. 
The sound speed in such a system is affected by several things. 
Firstly, the path length that a signal has to travel is not a straight 



where the £ is a displacement from the equilibrium position. This 
force can be then expressed as 




Fig. 6. Sketch of a sound wave propagation in a linear aggregate 
(a) and in a porous aggregate (b). 



line and thus with a given speed, longer distances result in lower 
effective sound speed. For a small RBD aggregate (approxi- 
mately 300 monomers), the effective distance was larger by only 
a factor of 1.5. Thus structure of an aggregate has a limited im- 
pact on the sound velocity by increasing the path length. The 
second factor that may have an effect is the tangential force. In 
a linear aggregate, the signal is passed forward due to the pres- 
ence of the vertical force, and the tangential force has no effect. 
However, in irregular aggregates grains interact also via rolling 
and sliding. This might lower the contribution of the vertical 
force and in this way change the sound speed. Fig. [6] shows a 
sketch of how a sound wave propagates in different aggregates. 
Indeed linear aggregates involve only stronger, vertical forces 
and thus the signal travels faster. Irregular particles also make 
use of the tangential, rolling force. 

In order to estimate the sound sp eed due to the rolling fric- 
tion, we use again the spring theory. Domi nik & Tielensl d 1995b 
gives the recipe for the rolling friction to be 



F ro u = kg, 



(32) 



with k = 6ny - a spring constant. We then apply this in eq. [28lto 
get 



c, = 1 1 - - )(2ro - S )J f" 7 , ■ 



(33) 



The dependence of the sound speed on a grain size, c s oc r Q 112 , is 
in this case much stronger than in the case of the vertical force 

/ -l/6\ 

(c s oc r ). 

When we apply eq. [33] to the 1 .2 /jm silica grains we get 
sound speed c s = 16.5 m/s. This is a factor of 30 lower than 
the velocity derived for the linear chain. This suggests that the 
sound velocity in a porous aggregate might be significantly dif- 
fer ent from the speed derived for a linear chain of grains. 

ISironol (120041) derived compressive and tensile strengths de- 
pendence on density, based on experimental data. That relations 
were later used to develop an SPH model of large, mm sized and 
larger aggregates collisions. The sound speed in an aggregate 
made of 0.1 pm ice monomers was calculated to be c s = y/E/p. 
He used values for the Young's modulus and the density to be 
E = 6 x 10 5 Pa and p = 0.1 g/cm 3 , respectively. The resulting 
sound speed is then c s = 77.5 m/s. 

When we apply our formula for the sound speed to 100 
aligned 0.1 pm ice grains, we get sound speed c s = 885 m/s. 
This is almost an order of magnitude higher. We have to keep in 
mind however that the speed in a linear chain may be consider- 
ably different to the one in a real porous aggregate. In the rolling 
cas e, the theore tical sound speed is also high with c s = 250 m/s. 

iTeiserl (120071) performed a lab experiment, where he mea- 
sured the sound speed in a RBD aggregate with filling factor 
= 0.15 and 1 .5 //m silica monomers. He hit the dust cake from 
below and measured the response at the surface with a force sen- 
sor. The measured velocity was c s = 30 + 4 m/s, very much con- 
sistent with our results when the signal is assumed to be trans- 
ferred through the rolling degree of freedom. 

In order to numerically derive the sound speed and further 
test our model we performed a simulation of an RBD dust cake. 
The cake was being pushed from one side and we determined the 
response time of different particles in the aggregate. We com- 
bined the positions of the particles with their response time, 
which results in an average sound speed in the dust cake. Fig|7h 
shows the result of our simulations. For each monomer in the 
dust cake we plotted a corresponding sound speed. Initially large 
spread in the data shows that inside the cake very linear struc- 
tures are present. This leads to a few particles with very large 
sound speed. At larger distances, however, the range of sound 
speeds is much lower and shows that a signal is transported 
mainly via rolling degree of freedom. The average sound speed 
that is determined at the far end of the dust cake is only a factor 
of about 1 .5 larger than the sound velocity in the rolling degree 
of freedom. 

We applied two different forces to the cake and later deter- 
mined mean sound speeds in the cake for both cases. When a 
stronger force of 5 x 10~ 6 N was applied to the cake the sound 
speed reached c s « 60 m/s, while with the lower force of 5 x 10 -8 
N the sound velocity was c s « 20 m/s. For one monomer the lim- 
iting sound velocity, as calculated for linear chain of monomers 
was overcome. This, however, happened in the case of larger 



■Froii = fwryf, 



(31) 
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Fig. 7. Sound velocity in m/s versus distance from the plane that 
hits the aggregate. Squares indicate applied force of 5 10" 8 N 
and triangles 5 10~ 6 N. The left panel (a) shows results of sound 
speed determination in a RBD aggregate, while the right one 
(b) in a CCP aggregate. The inset image shows the CCP ag- 
gregate placed in between two planes. The lines indicate sound 
velocities determined theoretically using vertical force in a lin- 
ear chain of monomers (dashed line) and using rolling friction 
(solid line). The Dashed area limited by dotted lin es indicates th e 
sound speed range determined experimentally bv lTeiserl d2007l) 



force, and thus the low displacement approximation was most 
likely violated leading to larger sound speeds. 

The simulations show that the sound propagation is domi- 
nated by the rolling friction. In order to verify this we run another 
simulation to determine the sound speed in the aggregate with 
monomers arranged according to the cubic close packing (CCP) 
because this arrangement disables any rolling motion. The inset 
in fig. I7J? shows the aggregate. We applied a force on one side 
of the agglomerate and determined the response time of differ- 
ent particles. Data shown in fig. [TJd shows that indeed closely 
packed aggregates are characterized by much higher velocities 
than the porous ones. The monomers in each layer received the 
signal at a different time because the compressing planes were 
simulated by big grains and thus central particles were hit first 
and they forwarded the signal. Thus particles at the edges of the 
aggregate responded later than the ones in the center of each 
layer. The velocity c s a 600 m/s shows that indeed sound speed 
in a porous medium strongly depends on porosity because of 
the forces involved in the transport of the signal and the longer 
path for the signal to travel in more porous aggregates. Using 
a constant sound speed in SPH simulations is bound to lead to 
spectacularly wrong results. 

4. Conclusions 

We ha ve modified the origina l code by iDominik & Tielensl 
d 19971) and IDominik & Nuboldl d2002l) . We have shown that a 
numerical N-particle method for studying the properties of ag- 
gregates can be calibrated to experimental results by including 
the flattening of small asperities on the surface of the grains, and 
by using critical displacements for rolling of grains consistent 
with measured values. The new code reproduces the measured 
critical velocities for sticking very well. We would like to em- 
phasize that there is currently no proof for plastic deformation 



actually happening at the grain surfaces, but as a model it works 
very well. 

We went on to measure the compressive strength of the ag- 
gregates and compared the results with experiments. We get very 
similar results with a small offset in porosity which is most 
likely due to the relatively small aggregates used in our simu- 
lations. Finally, we computed the sound velocity in an adhesion- 
dominated porous material and showed that this leads to very 
interesting results. Three very different velocities play a role in 
such a medium, and the different velocities are separate by fac- 
tors of the order of 10. The fastest speed is the bulk material 
sound speed which only plays a role after a body has been molten 
and re-hardened. The second speed, a factor of ~10 lower is the 
speed at which a signal is transported in a longitudinal wave 
in a linear chain. This speed applies either in a perfectly linear 
chain, or in an aggregate that has been compacted sufficiently so 
that rolling of grains is no longer possible. Finally, the slowest 
speed is the one transmitted by rolling forces in a non-straight 
chain of grains. A small decrease stems from the longer path 
the sound has to take in a porous aggregate. By far the largest 
fraction stems from the weak forces in the rolling degree of free- 
dom. Experiments show that this is indeed the dominant speed 
in porous aggregates. However, our results show that the sound 
speed should be a very steep function of density once a signif- 
icant number of monomers has 3 or more contacts with their 
neighbors. It is to be expected that SP H approaches to model the 
properties of dust aggregates dSironoll2004t ISchafer et alj|2007h 
will fail strongly if these effects are not taken into account prop- 
erly. 

In summary, we conclude that we do now have a working 
model of dust aggregates that can be applied in parameter stud- 
ies. 
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